Bulk viscosity and the phase transition of the linear sigma model 
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In this work we deal with the critical behavior of the bulk viscosity in the linear sigma model 
(LcrM) as an example of a system which can be treated by using different techniques. Starting from 
the Boltzmann-Uehling-Uhlenbeck equation we compute the bulk viscosity over entropy density of 
the LcrM in the large- N limit. We search for a possible maximum of £/s at the critical temperature 
of the chiral phase transition. The information about this critical temperature, as well as the 
effective masses, is obtained from the effective potential. We find that the expected maximum (as 
a measure of the conformality loss) is absent in the large-iV in agreement with other models in 
the same limit. However, this maximum appears when, instead of the large-iV limit, the Hartree 
approximation within the Cornwall- Jackiw-Tomboulis formalism is used. Nevertheless, this last 
approach to the LcrM does not give rise to the Goldstone theorem and also predicts a first-order 
phase transition instead of the expected second-order one. Therefore both, the large- N limit and the 
Hartree approximations, should be considered relevant and informative for the study of the critical 
behavior of the bulk viscosity in the LcrM. 
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I. INTRODUCTION 

Transport properties are essential to understand the equilibration of the hadronic matter created at heavy-ion 
colliders. In this context the most studied transport coefficient has been the shear viscosity rj. When two relativistic 
nuclei collide, the initial spatial anisotropy of the interacting region is converted into momentum anisotropy due to the 
£N| ', hydrodynamical evolution. The momcntum-anisotropy equilibration is mainly controlled by the shear viscosity and 
it influences the value of the flow coefficients v„ . Comparing the measured value of the elliptic flow v-i [H, with the 
numerical hydrodynamic simulations, a value of r//s ^ 0.1 has been estimated (where s is the entropy density). 

This small value, very close to the Kovtun, Son and Starinets (KSS) minimum conjectured in Q, has characterized 
the matter created in relativistic heavy-ion collisions as a nearly perfect fluid. It has also been shown that the KSS 
coefficient r//s has a minimum at the liquid-gas phase transition of common fluids @, A minimum near the phase 
transition of the linear sigma model (LcrM) in the large- TV limit has also been found @ and this minimum is also 
expected in the deconfinement phase transition of QCD jH. [Toj. 

Another transport coefficient that relates momentum flux with a velocity gradient is the bulk viscosity £. It is 
sensitive to uniform expansions of the system in such a way that it is closely related to the scale invariance of the 
fluid [Hldli]- For conformal systems the bulk viscosity identically vanishes. This coefficient has been assumed to be 
much smaller than the shear viscosity as it is the case for common fluids. Moreover, in perturbative QCD the ratio 
between the bulk viscosity and the shear viscosity has been estimated around 1CT 3 - 10~ 8 [13]. In the low temperature 
^3_" phase, the hadronic medium, this coefficient has been calculated using the Green-Kubo formalism and kinetic theory 



in [l4j and [l5(, respectively. In this phase, the ratio C/77 has been found to be around 1CT 3 - 3 • 1CT 3 [ltj . 

However, near the critical point C can be larger than the shear viscosity. For those systems belonging to the 
dynamical universality class of model H (like the liquid-gas transition) the bulk viscosity diverges [lj| with the 
correlation length as C - £ 2 ' 8 (or as a function of the reduced temperature [t = T/T c - 1] as C ^ t ). In lattice QCD 
calculations of the bulk viscosity over entropy density above the critical temperature [la , 1 1 91 ] this coefficient seems to 
diverge at T c . Other authors have pointed out that the bulk viscosity presents a maximum in the critical temperature 
for some different models like the O(N) model, even at mean-field level [20Tj23| . However, this maximum at T c is not 
seen in other systems as in the Gross-Neveu model in the limit of a large number of fermion fields [24| . where the 
bulk viscosity is a monotonically decreasing function when temperature increases. In the context of the O(N) model 
it has been recently shown [25[ that it belongs to model C if N = 1 and to model G if N > 1 (particularly the large- iV 
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limit) with a divergence with the correlation length if N = 1 going like £ ~ £ 2 and being finite for N > 1 as ( ~ £ . 

In the hydrodynamic calculations of relativistic heavy-ion collisions, the bulk viscosity over entropy density can 
also be extracted. For example in Q it is found that a value £/s = 0.04 is compatible with the integrated elliptic 
flow measured by RHIC and ALICE collaborations. In [26| a similar value is found (£/s < 0.05) near the freeze-out 
temperature. 

In principle our goal is to study the behavior of the bulk viscosity over entropy density in the LcrM in the large- N 
limit following the same lines of our previous article Q. We will compare the phase diagram of the LcrM in the 
large- iV limit with the coefficient £/s in order to ascertain if a maximum of £/s is found near T c . In Sec. [IT] we present 
the LcrM and its effective potential at finite temperature in the large- N limit. In Sec. Mil we review the calculation 
of the shear viscosity over entropy density in this model. In Sec. IIVI we will perform the detailed computation of the 
bulk viscosity over entropy density and discuss the results in terms of the conformality loss of the system. In Sec. |V] 
we repeat the calculation for £/s in the context of the Cornwall- Jackiw-Tomboulis (CJT) formalism in the Hartree 
approximation. Finally, in Sec. I VII we present our main conclusions. 



II. LcrM LAGRANGIAN IN THE LARGE- iV LIMIT 



In this section we will review the dynamics of the LcrM at finite temperature, deriving the effective potential in the 
large- N limit. As we have already detailed this procedure in Q we will only sketch the key steps of the calculation. 
The bare Euclidean Lagrangian of the LcrM reads [4!| 

c = ^d^ T d^ - tz 2 $ t $ + A ($ T $) 2 - e ** +1 , (i) 

where the multiplet $ = (iTi,cr) (i = 1, N) contains N+ 1 scalar fields. A is positive in order to have a potential bounded 
from below and we consider ~p 2 to be positive in order to provide a spontaneous symmetry breaking (SSB). When 
e = 0, the SSB pattern is SO(N + 1) -> SO(N). In the following we will refer for historical reasons to the tt fields as 
pions (even if it is well established that the LcrM is not an accurate model for chiral dynamics) and to the a related 
degree of freedom as the Higgs. 

The factor e = m 2 /^ is responsible for the physical pion mass and, as it is well known, it gives rise to an explicit 
breaking of the SO(N + 1) symmetry. In the limit T = 0, the potential has a nonzero vacuum expectation value (VEV) 
and one expects SSB. Choosing the VEV, which will be denoted by f w , pointing in the N + 1 direction, we get the 
equation: 

-2fU + ^fl-e = Q. (2) 

For small e the solution to this equation is: 



/- = \/^ + l=2 =A( £ = 0) + 1=2 • ( 3 ) 



'^ + -V/*(e = 0) + — a 
2A 47Z 2 v ' Ajl 2 

The VEV can be written in terms of the TV-independent F parameter defined as: 

(<f> T <f>) = (a 2 (T = Q)) = f 2 = NF 2 . (4) 

In our notation we will call f n (e = 0) and /„• the VEV at T = for the cases without and with explicit symmetry 
breaking term, respectively. In the next sections, the VEV at arbitrary temperature will be denoted by v(T), in such 
a way that = v(T = 0). 

At T = the low-energy dynamics is controlled by the broken phase. In this case (ir a ) = 0. Then the relevant 
degrees of freedom are the pions which correspond to the (pseudo) Goldstone bosons when (e + 0) e = 0. Fluctuations 
along the a direction will be denoted by a and they correspond to the Higgs, the massive mode which is relevant at 
higher energies (or temperatures): 

CT = fir + O -> (cr) = .U ■ (5) 

The Lagrangian ([1]) written in terms of 7r a , f n and a reads: 
C = -d^dV + -d^a - /tVtt - JI 2 d 2 + — [2f 2 (ir a ir a + 3a 2 ) + (ttV) 2 + 27r a 7r a c> 2 + a 4 + 4/^a 3 + 47r°7r%er] 



+ 



(-e + *±f* - 2m 2 /.) a - JI 2 f 2 + A/4 _ e/ „ (6) 
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Note that the tadpole term vanishes because of Eq. ([2]). The tree-level Higgs mass can be read from the Lagrangian 
in @ 

Ml = -27J 2 + , (7) 



and it can also be written as: 



M| = 4^ + 3f = M + _l,^(e = 0) e 

fir N f n N U K ' 



using Eq. ^ for the last identity. 
The pion mass reads: 



ml = -2-f + Qf* , (9) 
which, as expected, depends only on the explicit symmetry breaking term 

ml = f , (10) 

and vanishes when e = (obviously all pions are degenerate) . The pion and Higgs masses are related through 

Ml-ml + ^. (11) 
For small enough explicit symmetry breaking we have: 

T? - «^ , (12) 
A ■ {^V-'S-O ■ (13) 



Equation (|T3|) can be expressed in terms of the VEV at e = 0: 

A . » at*, (14) 

8fi(e = 0) a 2 
where a is a multiplicative constant relating and / w (e = 0) 

U=aU(e = 0), (15) 

which can be written as: 

Ml - 1ml (m s 
ot = — ^ — , (lo) 

where we have used Eq. Q and consider small e. 

In order to obtain the effective potential we will integrate out the fluctuations from the partition function Z. To 
do that it is convenient to introduce an auxiliary scalar field. This field allows for a systematic counting of TV factors 
and provides some simplifications in the large- N limit. Then, the integration of the pions and the Higgs is performed 
by standard Gaussian integration. 

The integration of the fluctuations is done regardless of their wavelenghts. All the frequency modes of the scalar 
fields are treated at the same footing and this "unorganized" integration produces two undesirable features in the 
effective potential. First, an imaginary part of the effective potential appears. This imaginary part has been given 
the interpretation of a decay rate per unit volume of the unstable vacuum state by Weinberg and Wu in (27j . 

The second characteristic is the nonconvexity of the quantum effective potential. However, the effective potential, 
defined through a Legendre transformation of the generating functional of the connected diagrams [28| , should always 
be convex. This nonconvexity problem and the imaginary part appear as long as a perturbative method is used to 
calculate the effective potential [2{|. 
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A possible solution for these problems can be given by using a nonpcrturbative method to generate the effective 
potential. For example, the functional renormalization group generates an effective potential in such a way that an 
organized integration of the fluctuations is performed. Following the ideas of the renormalization group, only the low- 
wavelength components of the quantum fluctuations are integrated out at each step. Thus, the UV components are 
infinitcsimally integrated step by step and the final effective potential (defined in the infrared scale) does not acquire 
an imaginary part and it remains convex at every scale (at the IR point, the Maxwell construction is dynamically 
generated by the renormalization flow) (30j . 

In spite of the previous discussion, we understand that it is not necessary to perform a more sophisticated method to 
obtain the effective potential. The only relevant information for us is the location of the effective potential minimum, 
which eventually gives the position of the critical temperature. This minimum appears always outside of the nonconvex 
region. On the other hand, the possible presence of an imaginary part (whose domain in fact coincides with the domain 
of the nonconvex part of the effective potential) is not relevant for us in this work. 



A. Auxiliary field method and effective potential at T + 

To compute the effective potential in the large- TV limit we start by considering the partition function: 

Z = J VTr a V<7 f d 4 x Cj , (17) 

with the Lagrangian of Eq. ([1]) . Then, we introduce an auxiliary field x m order to deal with the quartic coupling by 
using the Gaussian integral: 



(18) 



To prove the equivalence between the two Lagrangians note that this auxiliary field has introduced a mass term and 
a coupling with <f> 2 in the Lagrangian. However, it has no kinetic term which means that x nas n °t a t rue dynamics. 
The Euler-Lagrange equation for x simply gives x = 2\/2A<f> 2 /iV. Introducing this solution into the right-hand side of 
Eq. (fT5)l one obtains the Lagrangian 



c = -^ x 2 + — x$ 2 = --$ 4 + 2-9* = 

8A 2 N N N 



(19) 



which is the original interaction Lagrangian. 
The partition function can then be written as: 



Z = f V7T a VaV X exp - f d 4 x l-d^dV + -d^a - Ji 2 ir a ir a -Ji 2 a 2 - — X 2 + —X^^ + — X^ ~ ™ ) 
J J y 2 2 8A 2 2/ 



Thus the action in terms of the 7r a , a and \ fields reads: 

1 n { — 9 rrz \ n 1 



S = f d 4 x (- u E -27I 2 + V2 X ) 7r a + iff (- u E -2ji 2 + V2 X ) a- ^; 



X -to 



(20) 



Before identifying properly the pion propagator one must get rid of the unphysical a tadpole. We have already seen 
that this term vanishes at T = 0. Now, we perform a shift of the a field a = v + a which also produces a change in the 
auxiliary field x = X + 2\/2-^-v 2 + 4\/2j^va allowing to cancel the tadpole term for a and the unphysical mass mixing 
term between a and %■ After some manipulations we obtain the action: 



S= fd 4 x K a (- aE+ G-\0,x])7r a + ^(-n E+ G- 1 [0,x] + S^v 2 )a + ^v 2 G- 1 [0,x] 



N _ 2 y/2 2 A 4 



8A' 



where we have introduced the function: 



G n 1 [q, X ]^q 2 -2jl 2 + V2, 



X 



XV + — v - ev 
2 N 

(21) 



(22) 



Now comes our approximation for the auxiliary field, x is n °t going to be integrated out, but treat it at mean-field 
level, so that it contains no fluctuations. In particular, we apply this simplification for the quadratic term in the 
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action. At mean field (note that we will keep the same notation Xi instead of using (x) as they coincide in this 
approximation) : 



X = x-2x/2A„2. (23) 

So the quadratic term reads 

N N o V2 2 A 4 

-8A* ="8A* + - V *~N V ■ (24) 

The last two factors are cancelled in the action (|20|) . We are going to trade x by G~ [0, x]- Using (|22|) . the quadratic 
terms is: 

We will neglect the last term in the action because it is a constant. Finally we use the relations 

NJ1 2 / 2 (e = 0) H NF 2 



4A 2 2a 2 2a 



(26) 



to introduce F = v(T = 0), i.e. the VEV at zero temperature. 
The final action becomes: 



Stn*,v,a,G?[0,x\] = fdtx (-□ £ +G; 1 [0,x])7r a + ia (- o E +G; 1 [0,x] + 8^Vj s 

(27) 



1 A/" F 1 ^ AT 

+-« a G; 1 [0,x] - -^-G^O.x] - ^(G^O,*]) 2 - £W 



Notice that G 7r 1 [g,x] is nothing but the inverse of the pion propagator in Fourier space. The inverse propagator of 
the a field is: 

G a 1 [q lX ]=G- 1 [ qiX ]+&±v 2 . (28) 

In order to generate the effective potential for v we now integrate out the fluctuations. By performing a standard 
Gaussian integration of the pions we get: 



/ 2Vexp(- / d 4 x ±**[-n E+ G?[0, X ]]« a ) = f d 4 xcxp(| f r^log^.x]) , 



(29) 



where go = 2itnT is the well-known Matsubara frequency appearing in finite temperature computations. The same 
procedure can be applied to integrate out a. To be able to perform the integrations we have assumed that the auxiliary 
field (or G~ [0,x]) is homogeneous, i.e. it does not depend on x. This assumption is also taken for the VEV of the a 
field and allows to obtain a simple representation in Fourier space, where different modes do not mix between them 
and the integration is straightforward. 
The effective potential (density) reads 

1 / NF 2 \ N TV _c 1 _f 

V eff (v, G- 1 [0, x]) = -\v 2 -—j G-J [0, x] ~ j^iG- 1 [0, X ]) 2 -ev + - log G, 1 [q, X ] + - log G, 1 [q, X ] , (30) 



with 



(27T) 



Looking at the N power counting of the different terms in the effective potential one finds that all of them behave 
as O(N) except the last one. Therefore the contribution of the Higgs to the effective potential is suppressed in the 
large- N limit by one power of N and it will be neglected in the following. Then the effective potential becomes: 



V eff (v,G- 1 [0, X ]) = ^v 2 -^^G-\0, X ]-^(G-\0, X ]) 2 -ev + ^ £logG?[q, X ] ■ 



(32) 
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T (MeV) T (MeV) 



FIG. 1: Second-order phase transition. Left panel: Order parameter or VEV for a and its susceptibility. Right panel: Higgs 
and pion effective masses. 



The last term needs to be regulated because it contains a divergence. This divergence can be absorbed by a proper 
renormalization of the quartic coupling (9|. Thus the renormalized (/i-indcpendent) effective potential finally reads: 



K//( v ,G- 1 [o,x]) = i( v 2 -^)G;Ho,x]-™-f3o(r,G;Ho,x])-^(G x 1 [o,x]) 2 

where the function go{T, G^ 1 [0, x]) is defined in App. |A"1 
The stationary conditions are given by: 



J L log v^tm 

(33) 



dV e 



efj 



0. 



eff 



G; 1 [o,x]=G; 1 [o,x] 



dv 



, 



(34) 



and they provide G^gfOjx] and the order parameter vq in this approximation. The effective mass of the pion is 
obtained as m 2 = G^-gfO, x] and the effective Higgs mass as M R = G'^fO^] + which has the same form as 

Eq. (I]). 

For the e = case there is no explicit symmetry-breaking term and we expect to have a second-order phase transition 
defined by the critical temperature T c = \/T2F = v/l2/iV f„ . To obtain the numerical results appearing in Fig. [1] we 
have used a vanishing pion mass at T = (with TV = 3), a Higgs mass of Mr = 500 MeV and v(T = 0) = 93 MeV. In 
the left panel we show the behavior of v(T) that follows the analytical solution v(T) = f v \/l - T 2 /T^. The numerical 
critical temperature coincides with the theoretical value of T c = 2f n = 186 MeV. We also show the susceptibility 
defined as minus the T-derivative of the order parameter. Its peak shows the position of the critical temperature. 
In the right panel we show the effective masses as a function of the temperature. The effective mass of the pions at 
T <T C must be always zero according to the Goldstone theorem (numerically it is fixed at 0.5 MeV at T = in order 
to avoid computational problems). At T c it starts growing with temperature in the symmetric phase. The effective 
mass of the Higgs follows the same pattern as the order parameter becoming zero at T c . For higher temperatures it 
increases until being degenerate with the effective pion mass as expected because of the restoration of the SO(N + 1) 
symmetry. 

For the e t case a crossover instead of a real phase transition is expected (this is the same situation as adding an 
external magnetic field to a ferromagnet). We fix m T = 138 MeV at zero temperature and the same values for Mr 
and f v taken for the e = case. The results are shown in Fig. [5J The left panel shows the order parameter v(T) 
that decreases with temperature and never becomes exactly zero. The crossover temperature can be defined as the 
position of the peak in the susceptibility that we plot in dotted line. In the right panel we show the pion thermal 
mass (solid line) and the Higgs mass (dotted line). They are degenerate for high temperatures. 



III. SHEAR VISCOSITY OVER ENTROPY DENSITY 

The shear viscosity for the LcrM was obtained in [9|. In the following we will briefly review the method used there 
with some minimal changes. The transport equation for the one-particle distribution function f p (t,x) is: 



(35) 
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FIG. 2: Same as Fig. Q]but in the crossover case, when the pion mass is fixed at T = to = 138 MeV. 



Taking into account only elastic collisions (we will discuss later the influence of inelastic terms) this equation reads 
dfpi ^ x) + -|.v/ p (t,x) = |- fdT 12 , 3p [/i/ 2 (l + / 3 )(l + / P )-/8/p(l + / l )(l + / 2 )] , (36) 

where 



i 3 jr. 

dr 12 , 3p = — EI , 9 X30 F I 7 '" (2^)^ (4) (fcl + ^ " fc 3 " P) • (37) 

At first order in the Chapman-Enskog expansion the distribution function is expressed as the local equilibrium 
distribution function plus a small correction: 

f P = n p + fW . (38) 
where n p is the local Bose-Einstein distribution function 

«!»(*. x) = gj^jpgjj) ■ (39) 

e T(t,x) _ i 

The Boltzmann-Uehling-Uhlcnbeck (BUU) equation is then linearized in f^p: 

On (t x) NE r \ f (1) f (1) f (1) f (1) 

dr^ + iV^ r ^ (1+ni)(1+na) h v h h 

at 2 J n 3 (l + n 3 ) n p {l + n p ) n\(\ + n\) n 2 {l + n 2 ) 

(40) 

The left-hand side depends only on the space-time derivatives of n p , which can be explicitly obtained using the Euler 
and continuity equations. 

The shear viscosity can be also expressed in terms of f p ■ 

W« aN flj^ f P™> (41) 

where Vij = diVj + djVi - |<5^ V • V is the shear gradient of the velocity field. 

The left-hand side of the BUU equation also carries the same gradient (neglecting the influence of other transport 
coefficients) 

dn (t x) 

Ep {)t + P ' Vn f^' x) = ^ + n pW V v ■ ( 42 ) 
In order to cancel out this factor from the BUU equation and also from (|4ip the function f^p is taken as 

= -n p (l + n p )/3 3 B(p) p l p>V l3 , (43) 
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with B(p) being an unknown function of p. Inserting Eq. (|43[) into (|41[) one gets: 



30 



It is convenient to write this formula in terms of the adimensional variables 

E p m^ 
m n 1 

Introducing the integration measure d^n{x\y) defined in App.jB] we find: 

Nm 6 f 

This B(x) function can be expanded in terms of the polynomial basis P n (x) defined in App. [B] 



(44) 



(45) 



(46) 



B(x) = £ b n P n (x) . 
By projecting the BUU equation into the space generated by P n one gets 

N 
n=0 

where the functions K l are also defined in App. [B] The collision integrals read 

^ I n ( 2 ^)32"g ro ^ (27r) ^ (4)(fci +k2 ~ h ~ p){1 + ni)(1 + n2)nanp 



(47) 



(48) 



nl 4m 2 ' 



PiPj 



Pi{p) + r-Pi(h) 2^Pi(tfi) T^Pii^) 



m% m 

h 3i k 3 3 Ic^lc 1 ! h 2i h 2 



Hi 



At the lowest order in the expansion (|4"T| the shear viscosity can be written as 



307r 2 T 3 C'' ' 



(49) 



(50) 



where the pion effective mass m„ is now taken from the gap equations of the effective potential and 



nn 



Nn 2 
Am 2 .'. 



h~ 2 I ft ( 2 t)*2E m l r l 2 ( 2 ^)^ (4) (fci + k2 - fc 3 ~P)(1 + + "2)n 3 n F 



kuk 



kukij 



p l p ) 



u3il,3j J.HjLl.7 h 2i h 2 3 



(51) 



The scattering amplitude in the large- ./V limit has been extensively described in Q which includes the tree-level 
7r - 7r amplitude and its resummation in pion loops (as they are of the same order in N) . The finite pion mass gives 
rise to new couplings to the pion scattering 31] and the corresponding amplitude should be added to the chiral one. 
In the large- N limit, the s-channel scattering is the leading one. As the amplitude is 0(1/N), the total cross section 
is 0(1/N), but the average cross section (which is included in the collision integral) is 0(1/ N 2 ). 

Note that, as ~ 0(1/N), the shear viscosity is 0(N 2 ). This is expected [HJ since the shear viscosity is 
proportional to the inverse coupling constant squared, and this is suppressed by one power of N. 

The obtained numerical results for N = 3, m^{T = 0) = and Higgs masses Mr = 0.2, 0.5 and 1.2 GeV are shown in 
Fig. |31 They are similar to that appearing in our Ref. 0] where we have used a slightly different parametrization for 

fp 1 ^. We have found a minimum of rj/s for the three cases, always greater than the KSS bound l/(47r). However the 
exact position of the minimum depends on the value of Mr. 
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FIG. 3: Viscosity over entropy density in the LcrM at large N for different values of Mr. 




FIG. 4: Comparison between the minimum of the viscosity over entropy density and the position of the critical temperature. 
The minimum is located close to it but slightly below. 



To check whether the minimum of rj/s corresponds to the location of the critical temperature one must compare 
the previous plot with the order parameter. We show this in Fig. [4] at different values of F(T = 0). In the left panel 
we show the normalized value of v(T). The position of the critical temperature (where the order parameter first 
vanishes) depends linearly on F. The same behavior is followed by the minimum viscosity over entropy density shown 
in the right panel. The position of this minimum is close to T c , but not exactly there (as shown in Fig. U) but slightly 
below. Notice that T c is Mfj-independent in the large- N approximation considered here. 



IV. BULK VISCOSITY OVER ENTROPY DENSITY 



In this section we will perform the calculation of the bulk viscosity in the LcrM in the large- TV approximation. We 
start by considering both elastic and inelastic scatterings. Therefore, we do not introduce a pion chemical potential 
in the calculation, as the pion number is not conserved in principle. This fact simplifies the thermodynamics with 
respect to our previous work [l5| since, instead of using the isochoric speed of sound v n and the compressibility k^ 1 , 
one only needs to consider the adiabatic speed of sound 




where P is the pressure and e the energy density. 

Additionally, as we will use a quasiparticlc description of the scalar fields, we must introduce a nonvanishing term 
dm^ldT that enters in the left-hand side of the BUU equation: 

P^n p (x)\ c = Pn p {l + n p ) f ^ - E 2 v 2 s + Tm^v 2 s ) V • V . (53) 
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Here it is useful to define a new T-dependent parameter m that includes the derivative of the thermal mass: 

,dml 



m 2 = m 2 - T 2 



dT 2 



Then the left-hand side reads 



p M 5 AI n p (x)|^ = /3n p (l + n p ) 



P 2 / 2 , ~ 2\ 

— -v s (p +m ) 



V-V . 



The first-order Chapman- Enskog correction to the distribution function f p = n p + fp 1 " > is 

fp 1 ^ = ~n p (l + n p )/3 A(p) V ■ V . 

The linearized BUU equation reads 



(54) 



(55) 



(56) 



n p (l + n p ) 



V 2/ 2 , -2\ 

— -v s {p +m ) 



NE. 



- J dr 12t3p (l + n 1 )(l + n 2 )n 3 n p [A(p)+A(k 3 )-A(k 1 )-A(k 2 )]+C r 



(57) 



where we represent by C e ; the elastic collision operator and by Ci n the part of the collision operator including inelastic 
scattering (that we do not explicitly detail here). As we have not fixed the pion chemical potential, this term should 
be present in order to allow particle number changing processes. 

The stress-energy tensor in presence of a thermal mass is written as [23|, [2(| [33| : 



N 



J (2tt)3£ p X P 



2 dm 2 \ (1) 
dT2 I Jv » 



(58) 



that reduces to the usual stress-energy tensor when the mass is T independent. 

This new term does not modify the expression for the bulk viscosity in terms of fp 1 ^ because the former only 
depends on the spatial components of the stress-energy tensor and in the local rest frame one has u l = 0. Therefore 
we have: 



A N r d 3 p . . .p 2 
C= — / t — — nJl + n p )A(p)— . 
T J (2tt) 3 E v py v! u ' 3 



(27T) 



(59) 



However, this new term in the stress-energy tensor does change the form of the Landau-Lifschitz condition [HI r 00 = 
(notice that this is the only condition, since the one fixing the particle density number out of equilibrium does not 
apply here). 



r 00 = N 



r d 3 p ( 

J (2tt) 3 E p 



p 2 + m 2 )f^=0 



(60) 



As usual, this condition can be used to add a vanishing contribution to the bulk viscosity, making much easier the 
comparison with the BUU equation left-hand side (|55|) . 



N C d 3 p 
C= — / 7 — dr — nJl + n p )A(p) 
S T J (2tt) 3 Ep py v! yi ' 



P 2 / 2 , -2\ 

— - v s [p + m ) 



(61) 



Next, we consider again the same adimensional variables we have used for the shear viscosity (x = Ep/m^, y = 
rriTr/T) and also an integration measure dp^c,. This integration measure is described in detail in App. [5] including its 
corresponding scalar product, the norm, the moments, the functions P and the polynomial basis. The bulk viscosity 
is then expressed as the following scalar product: 



Nm 4 

C=^;(A(x)\P 2 (x)) c . 



(62) 



The projected BUU equation is obtained by multiplying both sides of Eq. (|57p by 4 ^ m 1 4 g Pi(x)d 3 p and integrating 
over the three-momentum: 

N r r 

{P(x)\P 2 (x)) = J d 3 p J dT 12 ^p(l + n 1 )(l + n 2 )n 3 n p P l (x)[A(p)+A(p 3 )-A(p 1 )-A(p 2 )] + (C m ). (63) 
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There is one important remark concerning the solution of the linearized BUU. In this linearized equation only one 
zero mode is present. When A{p) is proportional to x the right-hand side is zero due to the energy-conservation law. 
However, due to the presence of the inelastic collision operator, the zero mode associated to the particle conservation 
(A(p) oc 1) is absent. 

In the left-hand side it is easy to check that whereas {P1JP2) = we have (-P0I-P2) + 0, and therefore consistent with 
the previous remark. The BUU equation is solvable in the entire Hilbert space of solutions [HO] ■ 

Due to arguments of final phase-space and suppression in the large- N limit (that will be detailed later) , we will not 
consider the inelastic processes in the collision integral and retain only the 2 -*■ 2 processes. This simplification causes 
an inconsistency in the BUU equation, as the projection of the BUU onto P$ gives different results on both sides of 
the equation. This inconsistency is avoided by solving the equation in the subspace perpendicular to the "accidental" 
zero mode P = 1. 

Therefore we expand the A(x) solution: 

A(x) = £ ) ■ (64) 

71 = 1 

After symmetrization of the collision integral the BUU equation finally reads 

00 

fePMaOlf = X>«C„* , (65) 

n=l 

with 

AT— 2 r 4 j3,L. 

Cnl = ^J il (27r) 32fl. m 2 ^) 4 ^! + k 2 - k 3 -p)(l + m)(l + n 2 )n,n p A[P ro (x)]A[fl(g)] , (66) 

where A[P„(x)] = P n (x) + P n (x 3 ) - P n ( Xl ) - P n {x 2 ). 

For I = 1 one gets the identity = which does not determine the coefficient 01. The first nontrivial case corresponds 
to I = 2, for which the solution of the BUU equation reads: 

\\P2{x)\? , , 

a 2 = " 2 ; >U , (67) 

By introducing it in the formula for the bulk viscosity one finds: 

Nmi „„ , Nmi II 



C 



2tt 2 T 



In the following we proceed to the discussion of the results for the bulk viscosity over entropy density. They will 
depend on the value of the pion mass at zero temperature, i.e. on having a second-order phase transition or a crossover. 
We will begin with the latter. 



A. Crossover 



We start with the explicit symmetry-breaking case, for which we set a nonzero pion mass at T = 0. In order to 
get numerical results a possible choice could be the physical pion mass m^iT = 0) = 138 McV. The effective mass 
depends on temperature in such a way that m^iT} > m^iT = 0) = 138 MeV. For this reason, and analogously to the 
physical pion gas in [l5| . the inelastic terms in the collision operator (1 -*■ 3 or 2 -*■ 4) are suppressed by the Boltzmann 
exponential factor in the final phase space e~ 2m "/ T (the inverses being too improbable to occur if the gas is dilute). 

However, the crossover case does not have a precise definition of the "critical" temperature. Here we will define 
it as the point where the minus derivative of the order parameter (the susceptibility) peaks. From the left panel of 
Fig. [2] we obtain a value of T cr = 261 McV. 

In Fig. [5] we show the squared speed of sound for the pions together with the result from a pion gas with a constant 
(temperature-independent) mass of = 138 McV. The difference between the two curves is attributed to the effective 
mass of the quasiparticles. It is important to remark here that both results correspond to the noninteracting gas, 
i.e. ideal gas. The introduction of interactions in the thermodynamic functions (for example, through the free energy 
obtained from the effective potential) would be inconsistent with the conception of the Chapman-Enskog expansion 
at first order. 
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Pion m=m(T) 
Pion m=138 MeV 
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FIG. 5: v% in the LerM at large N for the crossover case (solid line) compared with the speed of sound of a pion mass with a 
constant mass of 138 MeV. 
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FIG. 6: (1/3- v%f as a measure of the violation of the conformality of the system. The bulk viscosity qualitatively follows 
this factor. 



The speed of sound turns out to be a monotonic function of the temperature. It takes the conformal value v s = 1/3 
at T = 259 MeV, very close to the crossover temperature. By Eq. (|B20[) one can check that at this precise point: 



T 



(69) 



When v 3 = 1/3 the squared norm of the source P2 identically vanishes. Therefore, the bulk viscosity is zero because of 
Eq. (|68| . More generally, the bulk viscosity features the factor (| - ) 2 which measures the violation of conformality 



in the system [231 ]. Thus, when the speed of sound takes values far from the conformal one, the bulk viscosity will 
be non- negligible, whereas if v 2 s is near the conformal value [5l|, the bulk viscosity will be close to zero. We plot the 
factor {\-v 2 s ) 2 in Fig. El 

Another measure of the loss of conformality is the so-called interaction measure. In D = 3 + 1 dimensions it is defined 
as (6)t = e - 3P, where e is the energy density of the gas and P its pressure. Using the relativistic equation of state 
P + e = Ts it can be written as (6}t - Ts ~ 4P as a function of the pressure and the entropy density. The interaction 
measure is shown in the left panel of Fig. [7J The result for the bulk viscosity over entropy density is plotted in the 
right panel of the same figure. As stated before, the qualitative shape of the £/s follows the behavior of (1/3 - v 
shown in Fig. |BJ 

A maximum in the bulk viscosity over entropy density is not seen at the crossover temperature. However, a 
minimum is found for v 2 s = 1/3. Later, we will discuss the apparent contradiction with the claim of a maximum of £/s 
appearing at the phase transition. 

To end this section, we want to stress that this result is perfectly consistent (with due respect to the inevitable 
differences) with the result in [24[ for the Gross-Neveu model in the large- N limit. In that case, no peak is seen in 
C/s. More specifically, in the crossover case (which is the equivalent to this one) the speed of sound is monotonic and 
structureless like in our case. However, the speed of sound docs not cross the conformal value because if D = 1 + 1 then 
vs\cft = 1 and this is the value to which the speed of sound approaches asymptotically at large temperatures [52| . 



2 \2 

s) 
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T (MeV) T (MeV) 

FIG. 7: Left panel: Interaction measure for the crossover case as a signal of loss of conformality in the system. Right panel: 
f/s in the LcrM at large N for the crossover case. 




FIG. 8: Left panel: v% in the LcrM at large TV for the second-order phase transition case. Right panel: Factor (1/3 -v%) 2 that 
measures the violation of conformality. 



Following our considerations, the bulk viscosity in that model should be a monotonic decreasing function on the 
temperature, going to zero at asymptotically higher temperatures (where the conformal value of the speed of sound 
is reached at T -* oo). 

In order to give a consistency check of the previous result we will repeat the calculation in the chiral limit (at zero 
temperature) for pions. 



B. Second-order phase transition 



If the pion mass vanishes at T = 0, then a second-order phase transition is expected. In this case we can neglect the 
inelastic processes by a different argument. In the broken phase, as the pion mass vanishes, there is no Boltzmann 
suppression in the final phase space. However, in the large-iV limit the tree-level amplitude for inelastic processes 
(2 -> 4) carries one extra 1/JV factor. The average cross section (or the average scattering amplitude inside the collision 
operator) is order 1/N for elastic scattering, but it is order l/N 3 for the inelastic case and therefore it is suppressed 
in the large- iV limit. 

Before looking at (/s we show the speed of sound in Fig. [8j In the broken phase the pions are massless and the 
speed of sound takes the value of a conformal gas w| = 1/3 (note that the Higgs does not contribute because its effect 
is suppressed in the large- TV limit). When the critical temperature T c = 186 MeV is reached, the mass of pions starts 
to grow and the speed of sound separates from the conformal value. In the right panel we show the factor of loss 
of conformality (1/3 - v|) 2 which gives an idea of the behavior of the bulk viscosity. Of course, we expect the bulk 
viscosity to vanish in the broken phase and to grow for T > T c . 

In the left panel of Fig. |H] the interaction measure is plotted. In the broken phase it also vanishes as the gas is 
conformal. The £/s coefficient is also shown in Fig. [5] The bulk viscosity is zero in the broken phase but it starts to 
grow and takes a monotonic behavior when temperature is increased. 
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T (MeV) T (MeV) 

FIG. 9: Left panel: Interaction measure as a signature of loss of conformality of the system. Right panel: f/s in the L<tM at 
large N for the second-order phase transition. 




FIG. 10: Top panels: Pion and Higgs thermal masses. Bottom panel: Order parameter. All of them are plotted as a function 
of the temperature and the pion mass at T = 0. 



Note that again this result is consistent with the Gross-Neveu model at large N in the chiral limit of fermion 
mass (24l |. In the high temperature phase the mass of the fermion field is exactly zero, and the speed of sound takes 
the conformal value of v s = 1. Therefore, the bulk viscosity turns out to be zero. No maximum is seen at the phase 
transition temperature. When decreasing the temperature (the fermion thermal mass increases), the bulk viscosity 
has a finite discontinuity and increases as T -*■ 0. This is exactly what happens in our case, but with the high and 
low-temperature limits reversed. 

Finally, we show how the second-order case connects with the crossover by slowly increasing the pion mass from 
0.5 MeV (our chiral value) up to 138 MeV. In Fig. [TU] we plot the evolution of the pion and Higgs masses and the 
order parameter v(T) as a function of the temperature and the pion mass at T = 0. In Fig. [IT] we show the squared 
speed of sound and the bulk viscosity. 
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FIG. 11: Left panel: Speed of sound as a function of temperature and pion mass at zero temperature. Right panel: Bulk 
viscosity over entropy density. Note the nonstandard direction of the axes for the sake of clarity. 
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FIG. 12: Squared speed of sound of pions with physical mass of 138 MeV. We use thermal chiral perturbation theory [39(] at 
different approximations. 



C. Discussion of the results 



In both cases, the second-order phase transition and the crossover, we have not obtained a maximum for the bulk 
viscosity at the phase transition. We have extensively studied the reasons for this fact. One of the keys to understand 
this result is the behavior of the speed of sound (or more generally, the equation of state). Conformality, or the 
lack of it, should be reflected in this factor cither by taking the conformal value v 2 s = 1/3 or going far from it. One 
can conclude that, in order to see a maximum of the bulk viscosity at the critical temperature, one should have a 
nonmonotonic behavior of v 2 s near T c . For instance, the bulk viscosity could approach the conformal value as T -> T c —, 
showing a sudden dip at T c and increasing again to the conformal value 1/3 for T >T C . This behavior would produce 
a maximum in the bulk viscosity. Actually, this is the scenario when the bulk viscosity is phcnomcnogically included 
in the QCD phase transition (see for instance [22|, [35| ) and from the lattice QCD equation of state [36l - l38j . 

This behavior is also consistent with the second peak in the bulk viscosity showed in [l4{ . The nature of this peak 
can easily be understood from the speed-of-sound curve. In that work a dilute gas of pions with a physical mass of 
138 MeV is considered. The speed of sound is obtained from finite temperature chiral perturbation theory at two 
loops and order T 8 . In Fig. [T^] we plot (dashed line) the results obtained by direct use of the formulas in [39| for the 
pressure of a pion gas. The local minimum of v 2 s is responsible for the breaking of conformality around T ~ 220 MeV 
and therefore a maximum in the bulk viscosity is found in [l4j at that temperature. 

Focusing again in the LerM in the large-iV limit, the monotonic behavior of the speed of sound is inherited by 
the behavior of the thermal mass of the pion. The mild increase of the pion mass, without any indication of phase 
transition, determines the dependence on temperature of the speed of sound and eventually the bulk viscosity. This 
behavior is a consequence of the large- N approximation, in which some of the interesting details of the model are 
washed out. In other words, the simplification of the large- N limit in some aspects of the model is also accompanied 
by an oversimplification of the dynamics that eventually gives rise to the absence of a maximum in the bulk viscosity. 
In fact, as we have already mentioned, from the dynamical renormalization group is not expected to have a divergence 
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of the bulk viscosity in the large- N limit of this model |25j . This limit was also considered for the Gross- Neveu model 
in [24[ and the results found there are perfectly compatible with ours. 

The use of a more complicated dispersion relation or a different thermal effective mass behavior can notably change 
this result. A hint for this fact can be found in [23|. In that reference the authors consider also the LctM model but 
in a different way. First, they do not take the large- N limit and therefore the interactions among pions and Higgs 
must be included. However, in order to avoid some divergences in the cross section, they only consider constant 
scattering amplitudes (in fact they do not take into account any correction to the amplitudes due to the finite pion 
mass). Second, they do not solve the BUU equation by the Chapman-Enskog expansion. Instead, they use the 
relaxation time approximation to obtain the bulk viscosity. Finally, the effective potential is calculated in the CJT 
approximation which gives qualitative different results on the behavior of the effective masses and eventually on the 
speed of sound and bulk viscosity. In that reference, a peak in the bulk viscosity is found in the crossover case for 
some specific values of the Higgs masses. 

Therefore, in order to conclude our study of the bulk viscosity in the LctM, it seems to be interesting to consider a 
different effective mass to show that a maximum in the bulk viscosity is obtained without changing any other aspect of 
our calculation. Thus we will use the Hartree approximation of the CJT theory as presented in [4fJ for the calculation 
of the order parameter and thermal masses. Doing so is not completely consistent since we are using the Hartree 
approximation for the effective potential and keeping the large- N limit in the scattering amplitudes. However, in the 
next and last section we will not try to develop a perfect consistent computation but just to introduce a different 
pion effective mass in our previous calculation to check if a maximum is obtained for the bulk viscosity. Additionally, 
we have also checked that the use of the scattering amplitudes in the large- N limit or the use of constant scattering 
amplitudes without any further correction give rise to similar results. 



V. THE CORNWALL-JACKIW-TOMBOULIS FORMALISM IN THE HARTREE APPROXIMATION 

In this section we introduce a thermal pion mass from the effective potential obtained within the CJT formalism for 
the LcrM in the Hartree approximation. Here we will refer to floj where this approach is nicely presented. In a nutshell, 
the CJT method is a tool for the computation of an effective action, not only for the VEV of the field 4>(x) = (Q(x)) 
but also for the two-point function G(x,y) = (T$>(x)$(y)). This generalized effective action can be understood as the 
generating functional of the two-particle irreducible graphs (2PI). Some of the two-particle irreducible graphs that 
contribute to the CJT effective potential are shown in Fig. [T5] 

In (4(j| two different approaches were used to sum up certain sets of diagrams, namely the largc-iV limit and 
the Hartree approximation. Of course the large- N approximation gives very close results to ours and it gives no 
new information about the bulk viscosity. Therefore, we will consider the Hartree approximation in which only the 
"double bubble" diagram is taken into account, and which is equivalent to summing up all the "daisy" and "superdaisy" 
diagrams in the 1PI effective potential. The essential features of the method are sketched in the Appendix [Cl 

The effective potential is a function of the order parameter <p and two dressed propagators G a and G v , one for the 
Higgs and one for the pion, respectively. From the effective potential V, one can obtain the two gap equations 

and also minimize the effective potential with respect to the order parameter. If the dressed propagators are written 
as functions of some effective masses G" 1 = k 2 + M 2 , then the three equations give the following nonlinear system for 
4>, M a and Af T : 

M 2 a = -2JI 2 + ^4> 2 + 4XF(AU) + ^F(M a ) , 

Ml = -27J 2 + § 2 + mp± F{A . U ) + fF(M a ) , (71) 
= [-27i 2 + f 4> 2 + ^F{M a )+A\F{M^)](j ) -e , 

where F(M) is defined in Eq. (|C24[) . This system can be solved numerically in order to obtain the effective masses 
and the order parameter that enter into our bulk viscosity calculation. 

Of course, a renormalization program should be properly performed here. The technical aspects of the renormal- 
ization of the model can be found in (40| and in the references therein. However, in this work we will deal only with 
the temperature-dependent part of the effective potential. 
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FIG. 13: Left panel: Order parameter and susceptibility in the CJT formalism in the Hartree approximation with masses 
Mtt = 138 MeV and M a = 600 MeV at T = 0. Right panel: Thermal masses in the Hartree aproximation for the same set of 
parameters. 




T (MeV) T (MeV) 

FIG. 14: Left panel: Square speed of sound for the case M n (T = 0) = 138 MeV and M CT (T = 0) = 600 MeV. A local minimum 
is seen near the crossover temperature. Right panel: Bulk viscosity over entropy density. The first maximum appears at the 
local minimum of the speed of sound and the second maximum at the maximum of «|, these points correspond to the zones 
where the loss of conformality is larger. 



A. Crossover case 



Within this new formalism we start again by considering the crossover case with a Higgs mass of 600 MeV and 
a pion mass of 138 MeV at T = 0. In Fig. [T3] we present the results obtained by solving the nonlinear system in 
Eq. (|71[) . In the left panel we show the order parameter and minus its derivative with respect to the temperature. 
They look qualitatively very similar to the results found by using the large- N limit with a crossover temperature of 
T cr = 230 MeV. In the right panel we show the thermal masses, very similar also to those found from the large- N limit. 
However, one can see that the pion thermal mass shows a plateau at T ~ 220 MeV with a derivative approaching 
zero in that region. This behavior produces an effect which is clearly seen in the speed of sound (where only the pion 
effects are included jH] ) in the left panel of Fig. Q3J Its shape is drastically different from the one coming from the 
large- N limit computation. Now a new minimum is obtained at T = 220 MeV and it corresponds to a small loss of 
conformality around this region. The subsequent maximum is also a signature of the loss of conformality. Among 
them, the conformal value 1/3 is reached. Therefore, one expects a zero bulk viscosity surrounded by two maxima. 
This is exactly the case in which the zero of the bulk viscosity resembles our zero in the large- N limit, as we show in 
the right panel of the same figure. 

Now we can proceed to increase the Higgs mass. The solution of the nonlinear system in Eq. (|71|) is shown in Fig. 1151 
where a similar result to the previous case is obtained. The main difference is that the transition is more abrupt with 
a nice peak in the susceptibility. Now the curve shows a crossover temperature of T cr = 260 MeV. The effect on the 
pion mass is seen in the right panel of that figure showing a clear zone around the T cr where the pion mass possesses 
a negative derivative. The implications of this behavior become clearer in the speed-of-sound curve that we plot in 
the left panel of Fig. 1161 The previous minimum at T cr is much deeper showing a larger breaking of conformality. 
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FIG. 15: Same as Fig. [T3] but with M„(T = 0) = 900 MeV. The order parameter goes faster toward zero and the susceptibility 
clearly indicates the position of the crossover temperature. The pion mass has a zone of negative derivative near the crossover 
temperature. 
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FIG. 16: Same as Fig. [14] but with M a (T = 0) = 900 MeV. At T cr the square speed of sound possesses an abrupt minimum 
showing a large violation of conformality. The bulk viscosity shows a clear maximum at Tcr • 



This minimum produces a peak in the bulk viscosity at the crossover temperature. This is the maximum of C/s that 
one expects at the phase transition and the one that is lost in the large- N approximation. 



B. First-order phase transition in the chiral limit 

Finally, we take the chiral limit for pions at zero temperature. In the left panel of Fig. [17] we show how in the 
chiral limit the phase transition is of first order, with a clear discontinuity in the order parameter at T c = 190 MeV. 
One feature of the Hartree approximation, as opposed to the large- TV limit, is that it does not respect the Goldstone 
theorem. This can be seen in the right panel of Fig. 1171 where the pion mass in the low-temperature phase is different 
from zero. In addition, both masses show discontinuities at T c . This jump is also seen in the speed of sound in 
the left panel of Fig. [15] The bulk viscosity presents a clear peak in the phase transition temperature with a finite 
discontinuity inherited by the first-order nature of the transition in the Hartree approximation. 
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FIG. 17: Same as Fig. [T3]but with M^(T = 0) = MeV. The discontinuity of the order parameter reveals a first-order phase 
transition. In the right panel one can observe that the Goldstone theorem is not satisfied for the pion mass within the Hartree 
approximation. 
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FIG. 18: Same as Fig. 1141 but with M^(T = 0) = MeV. The speed of sound inherits the discontinuity of the first-order 
transition. A maximum in the bulk viscosity over entropy density is still found at the critical temperature. 



VI. SUMMARY 



In this work we have extensively investigated the loss of conformality and the bulk viscosity in the LcrM. In the 
large- TV limit, where the dynamics of the Higgs field is washed out, the phase diagram has been calculated in order 
to pin down the location of the critical temperature. The bulk viscosity has been calculated in this limit and we find 
that it vanishes at the conformal points but it is different from zero otherwise. However, a maximum of the bulk 
viscosity is not found in the large- N limit, consistent with other systems in the same limit as in (24|. Moreover, we 
find this result quite natural as recent dynamical renormalization group calculation shows that in this limit the bulk 
viscosity remains finite (it does not diverge) in the critical point [25| . 

However we have shown that a maximum of the bulk viscosity can be found if one considers the Hartree approx- 
imation in the CJT effective potential. Just by obtaining the effective pion mass from this approach in our kinetic 
computation of the bulk viscosity, we find that a maximum is obtained at the crossover temperature for some appro- 
priate region of the Higgs mass. Moreover, a clear maximum of £/s is obtained at the critical point when the pion 
mass is set to zero at T = 0. However, this approximation fails to reproduce the Goldstone theorem of the LcrM and 
gives rise to a first-order phase transition instead of a second-order one. Consequently, the large- TV and the Hartree 
approximations considered here seem to be relevant and informative to fully understand the behavior of the bulk 
viscosity of the LctM at the critical point. 
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Appendix A: go(T,M ), gi(T,M ) and their relation with the moments of n p 

The function go(T, M 2 ) and its derivative gi(T, M 2 ) are needed for the temperature-dependent part of the effective 
potential in Eq. (|33|) . The function go(T,Al 2 ) is defined as the following integral: 

g^M^^fJdx^-y 2 )^^ , (Al) 

where y = M/T. Note that in Eq. (f3"3")l M 2 corresponds to the function G^fOjx]- The notation M 2 is used to suggest 
that this factor will turn out to be the pion effective mass squared when the gap equation is solved. 

The derivative of this function with respect to M 2 , appearing in Eq. (jM]), defines the function gi(T,M): 

(A2) 

J dM 2 ' V ; 



and can also be written as: 



T 2 r°° x/r 2 - v 2 



.91 

'y 



In the limit y -*■ 0, i.e. the conformal limit as the pion mass goes to zero, the two functions are given by: 



T 4 7T 2 T 4 

go T,0) = _r(4)C(4) = — - , (A4) 
cir* 45 

91 (T,0) = ^r(2)C(2) = . (A5) 

27T Z 12 



These functions can be related to the moments of the Bosc-Einstein distribution function, n p [42j . These moments 
read: 

2"ia a -an (Tjm) = N J j^^-n, p-V'-p"" , (A6) 



with E p = \Jp 2 + to 2 . They can be expanded in a tensor basis in terms of some coefficients I n ,k depending on the 
temperature. These coefficients read: 

Nm n+2 r°° ,1 
Z nfc (T,m) = 7 ™ „ , f dxx n - 2k (x 2 -l) k+1 l 2 —— , (A7) 
■ v ' {2k + 1)!!2tt 2 Ji K ' e» x -l 

where x = E p /m and y = m/T. They satisfy the recursion relation 

Z«+2,fc+i(r,m) = 1 (l n +2,k ~ m 2 In t k) 7 (A8) 

"f" o 

and can be related to some thermodynamical quantities in equilibrium, e.g. P = I2.1, e = ^2, 07 etc. 

Performing a change of variables in (|A1[) and (|A3|) . one can express these two functions in terms of the I n ,k integrals. 
In particular: 

g (T,M 2 ) = ^X 2)1 (T,M) , (A9) 
9l (T,M 2 ) = 1 I ,o(T,M) . (A10) 



21 



Appendix B: K 1 , /' and their relation with the auxiliary moments of n v 

In this work we have used two different integration measures for the shear and bulk viscosities, respectively. These 
measures are naturally given by the form of the viscosities as an integration over the distribution function. In terms 
of the adimensional variables defined in Eq. they explicitly read: 



( e yx - \y 



d f i n (x ] y)=dx (x 2 -l) 5 ' 2 -^ -2, (Bl) 



d^(x;y)=dx (x 2 -!) 1 ' 2 -^ - (B2) 

(e yx - 1) 

where clearly x e f2 = [l,oo). The moments of these measures are related to the respective source functions of the 
BUU equation. They have been defined as 



K' 



J dfirj x l , (B3) 

p = J dfi ( x 1 . (B4) 
From these measures one can define the scalar products: 

(/ I g)v = f d^ n {x;y)f{x)g{x) , (B5) 

(f\9)t = f u d^(x;y)f(x)g(x) , (B6) 

and the corresponding norms in the usual way. In addition, one can define the corresponding orthogonal polynomial 
bases that expand the space generated by the x % . For the shear viscosity we use a monic orthogonal polynomial basis 
defined as 

P (x) = 1 (B7) 
I\(x) = x-^ (B8) 

2 K0K3-K1K2 , K\ - K,K 3 

rb(x) = x h s x h =: (B9) 

V ' K$-K Q K 2 Kl-K K 2 V ' 

However, for the bulk viscosity we follow the same convention except for the polynomial P2(x), which is fixed to 
be the source function in the left-hand side of Eq. (|57|) . 



Po(x) = 1 (BIO) 
Pi(x) = x (Bll) 

: G-*y-HP> (Bi2) 

Notice that now (Pi | P%) = 0, but (Po \ P2) is different from zero, so this basis is not completely orthogonal. Further 
information about the integration measures, scalar products and polynomial basis can be found with more detail (also 
for the thermal conductivity coefficient, not studied here) in |4lj . 

The two functions K l and P are particular cases of more general coefficients expanding the "auxiliary moments" 
of the Bose- Einstein distribution function n p (42[. The auxiliary moments read 

jaia2 - an (Tj m) = N r f p (1 + ) p ^... p a n (B13) 

J (2ir 
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In an appropriate tensor basis these moments can be expanded in terms of the coefficients J n ,k depending on the 
temperature. By using our adimensional variables, these coefficients can be written as: 

Nm n+2 r°° , pV x 

piTI)i5p/ (B14) 

Some of them can be related to thermodynamical quantities. For example: 

-jr = n; — = e + P; = s . (Bl5) 

The coefficients satisfy the recursion relation 

Jn+2M1 = 7TT~~^ {Jn + 2,k ~ TT? J n . k ) . (BI6) 

It is not difficult to see that K % and P are related with the k = 2 and k = coefficients, respectively 

Nm 6+i 

J M ,2 = -^tf, , (Bl7) 

Nm 2+i 

Jifi = -^2-^ ■ ( Bl8 ) 

Using the formula (|Bl5[) we can express the entropy density and the heat capacity in terms of the integrals P: 

Nm 5 > T x _ m ds Nm 5 / dm T 



6tt 2 T 2 



^ (Bl9) 



When there is no chemical potential, the temperature is the only independent variable, and the adiabatic speed of 
sound and the isochoric one coincide. Then, they can be written as 

v 2 s = ^ = \ T h ~ d Jr T . (B20) 

which provides an universal relation for the speed of sound for any system that can be described in terms of free 
quasiparticles with masses depending only on the temperature and, in particular, for a free particle system. 



Appendix C: The Cornwall- Jackiw-Tomboulis effective potential 



In this appendix we briefly review the basic ingredients of the effective potential calculation in the context of the 
CJT formalism and the Hartree approximation. For more details we refer the reader to flol |43| - |46| . To introduce 
this technique we start with the derivation of the effective potential in the simple case of standard A$ 4 theory at zero 
temperature and later we will extend the calculation at finite temperature for the LcrM in the Hartree approximation. 



1. A$ 4 theory at zero temperature 



The classical Euclidean action for the A$ 4 theory reads 

f 



5[$] = f 

J X 



m ~< -2 , X x4 



2 d^{x)d^{x) + — + -*(x) 



where we have defined 



(CI) 



(C2) 



In the context of the CJT method, the starting point is a generating functional of one and two point Green functions 
depending on a local source J(x) and a bilocal one K(x,y), defined as: 
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Z[J,K]= f exp|s[$] + f J(x)$(x) + - /" $(x)iT(a;,2/)$(y) 

J . ^as 2 J x,y 

The corresponding generating functional for the connected Green functions, W[J, iT] is defined as 



(C3) 



Z[J,K] 



The expectation value of the field <p(x) = (<&(x))j,k and the connected two-point function G(x,y) in presence of 
the sources can be obtained through the functional derivatives of W[J, K]. 



5W[J,K] 



SW[J,K] 1 

- ~l<P{x)ny) +G(x,y)\ 



SJ(x) ^ v " 5K(x,y) 2 
This two-point function should not be confused with the tree-level propagator D(x,y), whose inverse reads 

5 2 S[<S>] 



D- L (x,y) 



(5$(.t)(5$(j/) 



For the action given in (|C1|) this becomes 

D- l {x,y) =-{u x - m 2 - ^(x) 2> jS(x - y) 



(C4) 



(C5) 



(C6) 



The two-particle irreducible (2PI) effective action can be obtained as the double Legendre transformation of 
W[J,K]: 



6W[J,K] 



SJ(x) 



J{x)- f 

Jx.y 



5W[W K { x,y) 



5K(x,y) 



r[<t>,G]=W[J,K]- f 

J X 

which, by using, (|C4=|) can be written as: 

T[<P,G]=W[J,K]- [ <j>{x)J{x)-\ f K(x,y)<l>(x)<Ky)-l f G(x,y)K(y,x) 

Jx Z Jx,y Z Jx,y 



The stationary conditions for the 2PI effective action read 
5T[4>,G] 



8<j)(x) 



-J(x) - f K(x,y)4>(y); 
■Jy 



— - = --K(x,y) , 

oG{x,y) 2 



(C7) 



(C8) 



(C9) 



which leads to the VEV and the dressed propagator when the sources are set to zero. 
Following Cornwall- Jackiw-Tomboulis, the effective action can be written as [43[ 



T[<f>, G] = S[<f>] - \ Tr log G- 1 - i Tr (D^G - 1) + T 2 [0, G] , 



(CIO) 



where ^[0, G] is the sum of 2PI diagrams with G as internal propagators. Some of these diagrams can be seen in 
Fig. [19] and all of them contribute to the C JT effective action. In the following we will refer to the first diagram as 
the double bubble diagram. 

As usual, the effective potential (density) is obtained by assuming to be constant, so that one gets: 



V(<t>,G)=U{<t>) + \ f\ogG-\k) + \ f(D- 1 (k)G(k)-l) + V 2 [cl ) ,G], 

Z Jk Z Jk 



(Gil) 



with U((f>) = m 2 (j) 2 /2 + A</> 4 /4! being the classical potential. 

Now the equations of motion of the effective potential are obtained by minimizing V with respect to <p aud G: 



SVi&G) 



SV(^G) 
5G 



, 



. 



(C12) 
(C13) 
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FIG. 19: Two-particle irreducible diagrams which enter in IV The solid lines correspond to dressed propagators G. 




FIG. 20: Dyson- Schwinger equation for the dressed propagator G (represented by wider solid line) in terms of the bare 
propagator D (narrower solid line) and a self-energy insertion. 



Here the last equation provides the "mass gap equation" for Go in terms of (j>. Substituting this Go into Eq. (|C11[) , 
one gets an effective potential in terms of c6 only. Then, by minimizing this effective potential with respect to <j>, one 
obtains the order parameter (VEV) 4>q. 

In fact, Eq. (|C13|) is nothing but the Dyson-Schwinger equation for the propagator: 



G _1 = D' 1 + G) , 

where we have defined the self-energy as the functional derivative of the 2PI diagrams: 

5V 2 ((f>, G) 



E(0,G)e2 



5G 



(C14) 



(C15) 



On the other hand it is possible to show that each 2PI diagram with G as internal lines, corresponds to an infinite 
number of 1PI diagrams with the bare propagator D as internal lines. 

As the most important example for our work here we will consider the Hartree approximation for the calculation 
of V2- Basically, it takes into account only the "double bubble" diagram of Fig. fl9l (or equivalently diagrams 0(A)). 
Then for the A$ 4 theory we have: 



V 2 



A 



f G(fc,, 

J k 



and the gap equation (|C13[) for the effective potential reads: 



(C16) 



(C17) 



which is depicted in Fig. [2TJ1 

From this equation one realizes that the double bubble diagram with G as internal lines is equivalent to the full 
resummation of daisies and supcrdaisies diagrams in Fig. [5T] with bare propagator as internal lines. 



2. LcrM at finite temperature 



The extension of Eq. (|C13|) to the LcrM with N pions and a Higgs is straightforward. In the Hartree approximation 
the graphs contributing to V2 (and their numerical prefactors) are those of Fig. [22] Moreover, at finite temperature 
the fc-integration should be replaced by Matsubara summations. Defining 

//^„,k)=T£ f-^L/^k) , (C18) 



where uj n = 2itTn, the effective potential reads: 
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FIG. 21: In the Hartree approximation only the 2PI diagram in the left-hand side is considered. It is equivalent to all daisy 
and superdaisy diagrams with the propagator D as internal lines. We show some examples of those diagrams in the right-hand 
side. 




3 2N N(N + 2) 



FIG. 22: "Double bubble" diagrams contributing to the V2 term of the effective potential in the Hartree approximation. Solid 
lines represent pions and dashed lines the a. Their numerical prefactors in Vi are also shown, 



1 r N r 

V(4»,G) = U(4>) + -J logG ff 1 (fc) + - J logG-^k) 



1 r N r 

- J lD-^k^G^k) - 1) + - J [D- 1 (k,4>)G v (4>,k)~l) + V 2 (ct>,G„,G a ) 



where V2 is given by 



V 2 (^G a ,G^ 



(C19) 



4 



fG a (4>,k)] +N(N + 2)±[f fj G v (<t>,k) 

2N± [G a (<f>,k) f G n ((j),k) . 
N Jb Jb 



'P J/3 

The two gap equations are obtained by minimizing the effective potential with respect to G„ and G 

12A 



G-\4>) = D a 1 (q>)+4\f^GM) + — f p G M) > 
G- 1 ^) = D-\<b) + 4(N + 2)\fG^)+^fGM) 



(C20) 

(C21) 
(C22) 



To solve the system one can use the ansatz for the dressed propagator: Gj 1 = k 2 + M 2 , with Mi being the effective 
masses. Minimizing the effective potential (|C19j) with respect to </>, one obtains the third equation that closes the 
nonlinear system leading to the three parameters <p (eventually called v), M a and M n . Introducing the explicit form 
of the bare and dressed propagators, the effective potential finally reads: 



-Jl 2 4> 2 + ^4> A ~e<t> + Q{M a ) + NQ(M 7r ) 



^-M^F(M a )-^[-2TI 



■% 2 



M?jF(M v ) 



A 



+ 3-[F(M a )] 2 + (N + 2)X[F(M 1T )] 2 + 2XF(M a )F(M v ) 



N 



where 



r r d k r d i 

F<yM "> = J k 2 + M 2 = J (2tt) 3 2e7. + J J2H 



d 3 k 1 1 



(2tt) 3 2E k J (2tt) 3 E k e^ E " - 1 



(C23) 



(C24) 



26 



and 

Q(M) . 1 1 ♦ M») . / + r / ^ log[l - .-*] , (C25) 

with Eh = \/ k 2 + M 2 . As it is well known, both integrals have a T = (divergent) and a finite temperature-dependent 
part. The renormalization of the effective potential is discussed in [40( and references therein. In particular, there 
exist some difficulties in the renormalization under the Hartree approximation. See references |47| and [481 ] for 
more details (HJ. To our purposes we only take into account the finite thermal contributions F(M) -* Fp(M) and 
Q{M)^Q p {M). 
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